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Abstract 

The renormalized quasiparticle-RPA is reformulated for even-even nuclei us- 
ing restrictions imposed by the commutativity of the phonon creation op- 
erator with the total particle number operator. This new version, Fully- 
Renormalized QRPA (FR-QRPA), is free from the spurious low-energy solu- 
tions. Analytical proof is given that the Ikeda sum rule is fulfilled within the 
FR-QRPA. 
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The well-known Random Phase Approximation (RPA) has found many applications 
in different fields of physics, from chemistry and condensed matter to relativistic field 
theory M. In nuclear physics RPA has been extensively exploited to model properties of 
the excited nuclear states that allows to calculate intensities of various nuclear reactions, 
including the double beta decay (see reviews [|]]). 

The neutrinoless double beta decay (0^/?/3-decay) , which violates the total lepton number 
by two units, is a sensitive low-energy probe for new physics beyond the Standard Model ||,[S| . 
The observation of the Oz//3/3-decay would give unambiguous evidence that at least one of the 
neutrinos is a Majorana particle with non-zero mass |J. The current experimental upper 
limits on the 0^/9/?-decay half-life impose stringent constraints, e.g., on the parameters of 
Grand Unification and supersymmetric extensions of the Standard Model. 

It is important that the same nuclear structure methods can be applied to calculate 
intensities of 2^/3/3-decay, which is a second order process allowed within the Standard Model. 
In this case the quality of the RPA calculations can be directly checked by comparison with 
the corresponding experimental data available at the moment for a number of nuclei 
Most of the nuclei are open-shell ones and one uses the quasiparticle version of RPA (QRPA) 
which takes into account nucleon pairing correlations. 

The QRPA is able to reproduce the experimental data ||, but with the strength of a 
particle-particle interaction which is very close to the point where QRPA collapses. This 
well-known drawback of the RPA is due to overestimation of ground state correlations when 
one uses the quasiboson approximation (QBA) for the bifermionic operators and the Pauli 
exclusion principle (PEP) is violated. The proximity to the collapse makes the entire QRPA 
scheme unstable and calculation results unreliable. To cure this problem, the so-called renor- 
malized QRPA (RQRPA) was invented ||, which takes into account PEP in an approximate 
way. Namely, the expectation values of the quasiparticle number operators are consistently 
treated to be non- vanishing in the RQRPA vacuum. The RQRPA does not collapse for phys- 
ical values of the particle-particle interaction strength. It has been shown within schematic 
models that by including PEP in the QRPA good agreement with the exact solution of the 
many-body problem can be achieved even beyond the critical point of the standard QRPA 
(see, e.g. ]10[ and references therein). The RQRPA has been used in previous studies of the 
double beta decay [MFIPI . 

Nevertheless, it has been shown within the modern versions of RQRPA that the model- 



independent Ikeda sum rule (ISR) is violated [|12| -[T^|. The Ikeda sum rule states that the dif- 



ference between the total Gamow- Teller strengths in the (3~ and (3 + channels is 3(N—Z) ||T^ 
This drawback of the model is quite serious, because the derivation of the ISR is based just 
upon two principles: closure relation and average conservation of the particle number. Viola- 
tion of the ISR immediately means violation of one (or both) of these important principles. 
The drawback has not been cured neither in the self-consistent QRPA (T^j or in second- 
QRPA ||13|| . The restoration of the ISR claimed in the Ref. |16| actually works only for a 
special one-level model and can not be generalized for the realistic case as shown in Ref. |14 



At the moment, the second-QRPA is able to reduce the violation of the ISR to the level of 
few percents only when contributions from the three (and more) boson states are neglected 
in the boson expansion of the bifermionic operators |13| . 

It can be conjectured that the closure relation is not fulfilled within the standard for- 
mulation of the RQRPA. The vacuum of the RQRPA contains quasiparticles. The RQRPA 
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takes into account the fact that the quasiparticles block creation of the bifermionic QRPA 
bosons but does not consider the possibility, that the quasiparticles in the ground state take 
part in the transitions. This leads to violation of the unitarity and underestimates the ISR 
which should be fulfilled exactly. 

In view of these remarks, modification of the phonon operator in order to include so- 
called scattering terms (describing transitions of the quasiparticles) is unavoidable if one 
wants to restore the ISR within RQRPA. This has been understood for some time and 
attempts have been undertaken to treat the terms as additional independent constituents of 
the phonon operator |]17|| , in spirit of the thermal extensions of RPA [T8j for description of the 
particle-hole excitations in heated nuclei. That leads to increasing of the dimension of the 
QRPA equation system and the appearance of low-lying energy roots originating from the 
transitions of the quasiparticles. Although such a consideration has allowed to restore the 
ISR for finite temperatures jyj, the low-lying spectrum has no physical interpretation when 
the even-even nuclei are considered. For instance, the model would predict a number of 2 + 
excitations below the first quadrupole state, that have never been observed in experiments 



(see also the discussion in Ref. [0). Authors of a recent paper [^(J have succeeded in 
removing these spurious states by introducing new quasiparticles but they have again used 
the QBA in terms of the new quasiparticles and formulated the QRPA with respect to the 
HFB vacuum. 

The above-mentioned examples show that a correct formulation of the RQRPA is quite 
a delicate task. The goal of this Letter is to introduce a version of the RQRPA for even-even 
nuclei which is able to overcome mentioned difficulties, i.e. it fulfills the ISR and at the 
same time is free from the spurious solutions. 

Let us start from the general basis of the RPA, which assumes that an excited state of 
the nucleus in question, with the angular momentum J and the projection M, is created 
by applying the phonon-operator Qj M to the vacuum state |0 PPA ) of the initial, even-even, 
nucleus with the same nucleon number A: 

| JM) = Qi M |0+ PA ) with QjmK pa ) = 0. (1) 

At this point one usually introduces the following ansatz for the phonon-operator Q\ M : 

Qjm = E [X { ^,j)A\W, JM) - Y (tt/j) A(tt', JM)] . (2) 

tt' 

where A( TT / j), Yu^n denotes free variational amplitudes, which are calculated by solv- 
ing the RQRPA equations; A^(rr',JM) and A(tt',JM) (r and r' denote single-particle 
quantum numbers including isospin projection) are the two quasi-particle creation and an- 
nihilation operators: 



A^tt', JM) = [44,1 = £ C^ mT]rimTi a\ m a\, m ,. (3) 



JM 

m T ,m' T 



The quasiparticle creation and annihilation operators a+ mr and a rniT are defined by the 
Bogoliubov transformation 



Arm. / V ~V T U-, 



(4) 
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where c+ mT (c Tmr ) denotes the particle creation (annihilation) operator. The tilde symbol 
indicates the time-reversal operation, e.g. a TmT = (— l) jT- " v a r _ ?TiT . 

However it turns out, that the most appropriate way is to analyze the phonon structure 
in terms of the particle creation and annihilation operators (instead of the quasiparticle 
ones). Let us first mention two quite general conditions which should be kept in mind while 
constructing the phonon operator. Namely, because the phonon operator is assumed to 
create a nuclear state of A nucleons acting on the + ground state of the initial nucleus (also 
of A nucleons) , the operator has to 

1) have good angular momentum J and the projection M (it has been taken into account 
in Eqs. ©,(§); 



2) commute with the total particle number operator A = N + Z = 4m T c rm T - Neutral 

r,m T T 

excitations can be further distinguished from the charge-exchange ones using their 
commutativity with N — Z. 



4<V 



The first rule tells us that the simplest building blocks for the phonon operator are 
together with their hermitian conjugates. Applying the second rule 



JM' 



44< 



JM 



eads to a further restriction of the phonon operator structure: 

Qjm = E [x iT r',j)C\rr', JM) - y^G^rr', JM) 



(5) 



with &{tt',JM) 



4<V 



JM 



and C(tt', JM) 



\J-M 



C(tt', J - M). For the neutral 



excitations r and r' have the same isospin projections, whereas for the charge-exchange ones 
- opposite. Therefore, in terms of the original particle creation and annihilation operators Q 
consists of scattering terms only. Going now into the quasiparticle representation, one gets: 

C ] {tt\ JM) = u t v t ,A\tt', JM) + u t 'V t A(tt' , JM) 
+ u t u t <B\tt' ', JM) - v t v t ,B{tt\ JM) 

= u t v t/ A\tt', JM) + v T u T ~A (tt, JM). (6) 

with the following definitions: B^(tt', JM) = a\a T i ^ for the quasiparticle scattering 

term and A^ = A^ + {u T iv T > — u T v T B^j I (v^t — vl) for the bifermionic operator A^ . 
Now we are able to rewrite @ in the form which is similar to (0): 

Qjm = E [x^^A^rr', JM) - F (rrV) l (rr', JM)] , (7) 

tt' 

where X = u T v T >x — v T u T /y, Y = u T v T /y — v T u T /x. Thus, instead of A A, bifermionic 
operators A^A are now the basic terms in which the RQRPA should be formulated. The 
key point is that the latter automatically contain the quasiparticle scattering terms which, 
however, are not associated with any additional degrees of freedom. That means that there 
are no spurious low-lying solutions in the present theoretical scheme. 

From this point we can follow the usual way to formulate the RQRPA |9|] substituting 
everywhere A by A. The forward- and backward- going free variational amplitudes X and 
Y satisfy the equation: 
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A B 
B A 



X 
Y 



U \ X 
QRPA -U)\Y 



(8) 



where 



"4 — (®rpa 
B = -(0 



A. 



RPA\ 



A, 



H F , A 



'RPAI i 

Qrpa) j 



(9) 



Hp is the nuclear Hamiltonian and the renormalization matrix U TT > is 



U TT > = (0 



RPA\ 



\®rpa) — & Ta <) T > a 'T) T 



A{Tr',JM),A\aa',JM) 
V TT , = 1 + - v^/ )Af T > - (u 2 T - v 2 )Af T ) /(v 2 T , - v 2 T ) (10) 

Here, the exact (fermionic) expressions of the commutators are taken into account and 

K = j; 2 (0 + RPA \J2^m T arm T KpA) = ~Jr' Kp A \B(tT, 00)\0 + RpA ) (11) 

gives the relative quasiparticle occupation number for the level r in the RQRPA vacuum (j = 
\J1j + 1). The method to calculate these occupation numbers can be found elsewhere 0. 
It is useful to introduce the notation: 



X = U 1/2 X, Y = U 1/2 Y, 

A = U-^AU- 1 ' 2 , B = U-^BU- 1 ' 2 . 
Then the amplitudes X and Y satisfy the equation of usual QRPA: 



A B 
B A 



X 
Y 



1 

QRPA I g _j 



X 

Y 



(12) 



(13) 



(14) 



Solving the FR-QRPA equations, one gets the fully renormalized amplitudes X, Y with the 
usual normalization and closure relation, which will be used below to prove the fulfillment 
of the ISR: 

E™ -y-k \rm \rk r 

A (tt',J) A (tt',J) ~ I (tt',J) I (tt',J) ~ "km, 

TT 1 

(tt' , J) (t x , J) ~~ 1 (tt' ,jy (t^J) - Ott^t't'^ 



Evm \rra vm V m — n 

A (rr',J) r (r 1 r^,J) ~ (tt' ,J) (t^t^J) ~ U " 



(15) 



Here, m and k mark different roots of the QRPA equations for a given J 71 ". 

Now we are ready to prove the fulfillment of the ISR within the FR-QRPA. Let restrict 
ourselves to the charge-exchange pn-FR-QRPA and write down the definition for the ISR: 



ISR = £ \{JM,m\4 M KpA) ~ E \(J - M,m'\(3 JM K PA ) 



(16) 
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Here, we would like to consider the Fermi and Gamow- Teller transitions simultaneously by 
defining: 

Pjm= H (pm p \qjM\nm n )clcn = - J -1 ^(p||gj||n)C t (pn, JM), (17) 

pn,m p m n pn 

where qoo = 1 (Fermi transitions) and qi n = a u (Gamow- Teller transitions) and one has 
for both cases ISR = N — Z (in the latter case - just because of the only fixed angular 
momentum projection). One can use Eqs. (^),(0) to get: 

pn 



(J - M,m\f3 JM K PA ) = -(-) J - M J- ly E(n\\qj\\p) (u p v n Yr + v p u n XZ n , n Jv pn . (18) 



pn 



Substituting (|18|) into (|16|) and making use of the closure conditions (|l^), one ends up with 

ISR = J' 2 J2\( P \\qj\\n)\ 2 (v 2 n - v 2 p )V pn . (19) 



pn 

Within the usual QBA one has T> pn = 1 and 



ISR = J' 2 \(p\\qj\\n)\ 2 v 2 n - £ \(p\\<u\\n}\ 2 v 2 ) = N - Z, (20) 

\ np np / 



making use of J2 I (p||<?j|| n ) I = J Jn> J2\(p\\q.j\\ n )\ = J 31 an d of the equations for the 

p n 

BCS chemical potentials A n , X p 

Kfb\nK fb ) = YX< = N > Kfb\zK fb ) = £1V = z - (2i) 



n 



This corresponds to the fulfillment of the ISR for the usual QRPA. 

In the case of the FR-QRPA with V pn determined by Eq. (|T0"D one has 



ISR = J2j 2 n (v 2 n + (u 2 n - v 2 n )Af n ) 

n 

-HjI{v 2 p + {uI-vI)n p ) 
p 

= N - Z (22) 
because of the modified FR-QRPA equations for the chemical potentials: 
(0+ PA |iV|0+ PA ) = £ j 2 (v 2 n + (u 2 n - v 2 n )N n ) = N, 

n 

Kpa\Z\0 + rpa ) = E3 P 2 (*J + («J - v 2 p Wp) = Z. (23) 
p 

In conclusions, we have reformulated the RQRPA into Fully-Renormalized QRPA (FR- 
QRPA) for even-even nuclei using restrictions imposed by the commutativity of the phonon 
creation operator with the total particle number operator. The FR-QRPA is free from the 
spurious low-energy solutions. We have also shown analytically that the Ikeda sum rule is 
fulfilled within the FR-QRPA. 
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